Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 29 October 2012 (MN I^TeX style file v2.2) 



A Deep Chandra Observation of the AGN Outburst and 
Merger in Hickson Compact Group 62 

D. A. Rafferty/ L. Birzan,^ P. E. J. Nulsen,^ B. R. McNamara,^'^'^ 
W. N. Brandt,^ M. W. Wise*^ and H. J. A. Rottgering^ 

^Leiden Observatory, Leiden University, P.O. Box 9513, 2300 RA Leiden, The Netherlands 

Harvard- Smithsonian Center for Astrophysics, 60 Carden St., Cambridge, MA 02138, U.S.A. 
^Department of Physics and Astronomy, University of Waterloo, Waterloo, N2L 2C1 Ontario, Canada 

Perimeter Institute for Theoretical Physics, Waterloo, N2L 2YS Ontario, Canada 
^Department of Astronomy and Astrophysics, Pennsylvania State University, University Park, PA 16802, USA 
^ASTRON (Netherlands Institute for Radio Astronomy), P.O. Box 2, 7990 AA Dwingeloo, The Netherlands 



29 October 2012 



ABSTRACT 

We report on an analysis of new Chandra data of the galaxy group HCG 62, well 
known for possessing cavities in its intragroup medium (IGM) that were inflated by 
the radio lobes of its central active galactic nucleus (AGN). With the new data, a 
factor of three deeper than previous Chandra data, we re-examine the energetics of 
the cavities and determine new constraints on their contents. We confirm that the 
ratio of radiative to mechanical power of the AGN outburst that created the cavities 
is less than 10~^, among the lowest of any known cavity system, implying that the 
relativistic electrons in the lobes can supply only a tiny fraction of the pressure required 
to support the cavities. This finding implies additional pressure support in the lobes 
from heavy particles (e.g., protons) or thermal gas. Using spectral fits to emission in 
the cavities, we constrain any such volume-filling thermal gas to have a temperature 
kT ^4.3 keV. For the first time, we detect X-ray emission from the central AGN, 
with a luminosity of L2-10 keV = (1-1 ±0.4) x 10^^ erg s~^ and properties typical of a 
low-luminosity AGN. Lastly, we report evidence for a recent merger from the surface 
brightness, temperature, and metallicity structure of the IGM. 

Key words: X-rays: galaxies: clusters - galaxies: clusters: individual: HCG 62. 



1 INTRODUCTION 

In the past decade, images from the Chandra X-ray Obser- 
vatory of the hot atmospheres of clusters and groups have 
been critical in furthering our understanding of AGN feed- 
back. These images often show the radio lobes of the AGN 
displacing the hot gas, leaving cavities in the X-ray emission 
that are filled with radio emission (e.g., |McNamara et al.' 



200Q| [Fabian et al.|20QQ ; Blanton et al.|2001| ). These images 



are direct evidence of the strong coupling between the ra- 
dio lobes of the AGN and the hot, cooling gas. Estimates 
of the power required to produce the cavities has revealed 
that, in many such systems, the total jet power of the AGN 
is typically factors of 100 larger than the radiative power at 
monochromatic radio wavelengths (Birzan et al. 2004). In 



massive cooling flows, the total jet power is on average sim- 
ilar to the total cooling luminosity of the hot gas ( [Rafferty 



et al.|2006 ), whereas in lower mass groups and ellipticals, it 



tends to be factors of several larger than the cooling lumi- 
nosity, but with a lower duty cycle (e.g., Nulsen et al.|2009 



Dunn et al.||2010| [Pong et al.||2010| ). Therefore, it is criti- 
cal to study cavity systems in detail to understand fully the 
impact of the AGN on their group- and cluster-scale envi- 
ronments. In particular, the properties of AGN cavities at 
the group scale are less well known, due to the relatively few 
cases of clear AGN cavities in groups. 

We report on the results of an analysis of deep Chan- 
dra data on HCG 62, one of the best known cases of AGN 
feedback in a group. HCG 62 is a nearby {z = 0.0137), 
compact group of ^ 60 galaxies. Three large galaxies (NGC 
4761, NGC 4776, and NGC 4778) dominate the core, with 
NGC 4778, an SO galaxy, being the brightest cluster galaxy 
(BCG). The BCG hosts a low-lu minosity AGN with weak 
line emission ( Coziol et al.|[l998l ,Shimada et aLl[200Qt and 
weak radio lobes ( Gitti et al.|2010 Giacintucci et al.|2011 ). 
Mid-infrared Sptizer observations by Gallagher et al.|( 2008 ) 



show no evidence for excess nuclear mid-infrared emission, 
and no X-ray emission from the AGN was detected by 
orita et al. (2006), who found an upper limit for this emis- 
sion of Lo.5-4 keV ^ 10^^ erg s~^. Optical observations of the 
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BCG show no evidence of recent star formation (Rafferty 
et al.|2008| ). The group is X-ray luminous, with a 0.5-8 keV 
luminosity of ^ 10^^ ergs s~^. The temperature and den- 
sity profiles of the IGM (e.g., |Morita et al.|[2006t show the 
classical cooling-flow behavior in which the temperature de- 
creases inwardly to a cool, high-densit y core. T he cooling 
time in the core is ^6x10^ yr (R afferty et al.|2 006). 

HCG 62 was first observed by Chandra in 2000 for ^ 
50 ks. This observation revealed cavities in the IGM that 
were subsequently shown to be filled with radio emission 
from the lobes of the central A GN (fVrtilek et al.|2 002 ; Birzan 
et al. 120041 IMorita et al.|2006| |Gitti et al.|2010]]Giacintucci, 



et al. 112011 ). Interpreting these cavities as buoyant bubbles. 



Birzan et al.| (|2004|) calculated that the cavities represent a 

^, more than 



3 X 10^^ erg s" 



total mechanical power of 
enough to offset cooling losses of the IGM in the core. HCG 



62 has also been observed by XMM- Newton (Morita et al. 
2006||Gitti et al.|2010t, A SCA ( [Nakazawa et al.||2007t , and 
Suzaku ( Tokoi et al.|2008| , providing a wealth of X-ray data 
at a variety of scales, resolutions, and energies. Therefore, 
HCG 62 is an ideal candidate for detailed studies of various 
properties of the hot atmosphere in groups. 

We report here on new, deep Chandra observations of 
HCG 62 that reveal further cavities and evidence of a recent 
merger. We adopt Ho = 70 km s~^ Mpc~^, Qa — 0.7, and 

= 0.3, implying a luminosity distance to HCG 62 of 
59.3 Mpc and a angular scale of 0.28 kpc arcsec"^. 



2 DATA REDUCTION AND ANALYSIS 

HCG 62 was observed by the Chandra X-ray Observatory 
on 2000/01/25 for 48.5 ks (ObsID 921, FAINT mode), on 
2009/03/02 for 67.1 ks (ObsID 10462, VFAINT mode), and 
on 2009/03/03 for 51.4 ks (ObsID 10874, VFAINT mode), 
all with the ACIS-S detector. The data were reprocessed 
with CIAO 4. fusing CALDB 4.2.lQand were corrected for 
known time-dependent gain and charge transfer inefficiency 
problems. For each observation, blank-sky background files 
were used for background subtraction]^ The events files were 
filtered for flares using the CIAO script Ic-clean to match 
the filtering used during the construction of the background 
files. A total of 7.5 ks was removed during filtering, resulting 
in a combined exposure time of 159.5 ks. The background 
files were normalized to the count rate of the source image in 
the 10 — 12 keV band (after filtering). The derived normal- 
ization ratios (in the sense of background/source) were 1.01 
for ObsID 921 and 0.92 for both ObsID 10462 and ObsID 
10874. Lastly, point sources detected using the CIAO tool 
wavdetect were removed. 



2.1 Imaging Analysis 

Images used in the analysis that follows were made with 
the filtered events data. To correct for exposure, which is a 
function of both position and energy, images were made in 
12 narrow energy bands from 0.5 keV to 7 keV. The bands 



^ See cxc.harvard.edu/ciao/index.html 
See ^cxc . harvard . edu/caldb/ index . htinl| 



were chosen to enclose energy ranges for which the ACIS- 
S detector has approximately constant response. Exposure 
maps were weighted for each narrow-band image using a 
monochromatic energy equal to the central energy of each 
band. The narrow-band images and exposure maps from the 
three observations were reprojected to the same grid and 
summed together using the CIAO tool reprojectAmage. Re- 
gions where point sources were removed from the events data 
were filled with background emission using the CIAO tool 
dmfilth. The narrow-band images were then divided by their 
associated exposure map and summed together to produce 
the final, broad-band (0.5-3 keV), exposure-corrected image 
used in subsequent imaging analysis. 

The exposure-corrected 0.5-3 keV image is shown in 
Figurep^. A numb er of features noted by other authors (e.g., 
Morita et al.|2006 Gitti et al.|2010 ) in analyses of the orig- 
inal data are apparent, including the large cavities, a sharp 
inner surface brightness edge just beyond the cavities, and 
an outer surface brightness edge some 30 kpc (^ 2 arcmin) 
to the south of the core. 

We derived an azimuthally averaged surface brightness 
profile from the exposure-corrected 0.5-3 keV image by ex- 
tracting the mean surface brightness in 1-pixel-wide circu- 
lar annuli. No masking of the features mentioned above 
was done. The resulting profile is shown in Figu re [g] To 



parametrize the profile, we fit a double-/^ model ( [Sarazin 
BahcaII||1977 ) with a constant background to the profile 



The best-fitting parameters are given in Table [l] Although 
this model fits well over most of the profile, a clear depar- 
ture from the data is apparent at a radius of ^ 120 pixels 
60"), which corresponds to the location of the inner sur- 
face brightness edge. We find no evidence for a point-like 
source, indicative of X-ray emission from an AGN, at the 
centre of the profile, confirming the findings of [Morita et al.| 
( ,2006 ). 

We also fit a 2-dimensional circular double- /3 model to 
the core region only (r < 15 kpc). We allowed the positions, 
normalizations, and /3 parameters of the /3-model compo- 
nents to vary freely. Fitting was performed in Sherpa ver- 
sion 4.^ to the exposure-corrected 0.5-3 keV image, with 
the cavity regions masked. The best-fitting parameters are 
given in Table [l] These parameters differ somewhat from 
those determined from the large-scale azimuthally averaged 
surface brightness profile, but provide a better fit to the 
core emission. This 2-dimensional model is later used in the 



determination of the cavity sizes (see {3.2). 

To minimize artefacts due to asymmetries in the large- 
scale background group emission that are not accounted 
for in the /3-model fits, we also modeled the exposure- 
corrected 0.5-3 keV image using the Mult i- Gaussian Ex- 
pansion method of Cappellari et al. (2006). Although this 



method was developed for use with optical images of galax- 
ies, we have found it to work well with X-ray images that 
contain a large number of counts. We restricted the axial 
ratio of Gaussians to be between 0.2 and 1 to exclude highly 
elliptical Gaussians that are not appropriate for the smooth, 
large-scale background group emission. The residual image 
for this model highlights asymmetries in the distribution of 
the hot gas and is shown in Figure [TJ). 



^ See http: //asc .harvard. edu/contrib/maxim/acisbg/ 



See http : //cxc . harvard . edu/ sherpa/ index . html 
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Table 1. Parameters of the ^^-model Fits. 





Centre Position (J2000) 






Core Radius 


Comp. 


RA 


Dec 


Norm" 




(pixels) 


Azimuthally Symmetric Model 




12:53:05.7 


-09:12:14.1 


(1.24 ±0.01) X 10-'^ 


0.563 ±0.001 


35.3 ±0.1 




12:53:05.7 


-09:12:14.1 


(1.23 ±0.04) X 10-'^ 


2.23 ±0.07 


27.5 ±0.5 


BG 






(7.87 ±0.05) X 10-10 






2-D Model 




12:53:05.7 


-09:12:14.6 


(2.74 ±0.13) X 10-'^ 


0.33 ±0.05 


6.2 ±0.7 


/52 


12:53:06.0 


-09:12:09.7 


(9.0 ±0.65) X 10-8 


0.60 ±0.04 


44.2 ±2.5 



"Normalizations are for the 0.5-3 keV band and have units of counts s ^ cm ^ pixel ^. 
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Figure 2. The azimuthally averaged surface brightness profile 
(0.5-3 keV), with the best-fitting double- /3 model (solid lines) 
plus constant background model (dotted line) overplotted. The 
sum of the (3 model and background components is shown by the 
dashed line. The sharp inner edge is visible at a radius of ~ 120 
pixels (~ 1 arcmin). 



Lastly, to highlight the cavities and surface-brightness 
edges, we made unsharp-mask images at two scales: an image 
at large scales was made by subtracting an image smoothed 
at 32^^ from an image smoothed at 6^^ and an image at 
small scales was made by subtracting an image smoothed at 
from an image smoothed at 4" . These images are shown 
in Figures [l]: and[T]l. 

There are a number of features in these images that 
were not readily apparent in the original, shallower Chandra 
data. First, on large scales, one can see a spiral- like region of 
excess emission that extends to the south of the core, which 
has an outer boundary approximately coincident with the 
surface brightness edge in that region. The contours shown 
in Figure[T]3 highlight this excess. On smaller scales, in addi- 
tion to the large cavities found previously, two small cavities 
are now apparent on either side of the core, oriented roughly 
north-south. These features will be discussed further in ^ 
and g] 



2.2 Spectral Analysis 

For the spectral analysis, spectra were extracted from each 
observation separately using the CI AO script specextract. 
For each spectrum, weighted responses were made, and a 
background spectrum was extracted in the same region 
of the CCD from the associated blank-sky background 
file. Joint spectral fitting, using all three observations, was 
then performed when deriving spectral properties. For the 



spectral fitting, XSPEC |Arnaud]|T996| version 12.5.1 and 
Sherpa version 4.2 were used. 



2.2.1 Azimuthally Averaged Properties 

X-ray spectra were extracted in circular annuli centred on 
the centroid of the cluster emission. The annuli were con- 
structed from ObsID 921 such that each annulus contained 
~ 4500 counts, with the exception of the central annulus, 
which contained ^ 1000 counts (to allow derivation of the 
various properties as close to the core as possible). Spectra 
were then extracted from ObsIDs 10462 and 10874 from the 
same regions and generally contained a similar number of 
counts. Gas temperatures and densities were found by depro- 
jecting the se spectra using the Dire ct Spectral Deprojection 
method of Sanders & Fabian ( |2Q07[ ) Q The deprojection was 
done for each observation separately. The three deprojected 
spectra in each annulus (one for each of the observations) 
were then fit jointly in XSPEC with a single-temperature 
plasma model (MEKAL) absorbed by foreground absorp- 
tion model (WABS), between the energies of 0.5 keV and 
7.0 keV. In this fitting, the redshift was fixed to z = 0.0137, 
and the foreground hydrogen column density was fixed to 
A^H — 3.03 X 10^° cm~^, t he weighted- average Galactic value 
from [Dickey k Lockman] ( [T99Q| )p] 

The density was then calculated from the normalization 
of the MEKAL component, assuming Ue = 1.2nH (for a 
fully ionized gas with hydrogen and helium mass fractions 
of A = 0.7 and Y = 0.28). The pressure in each annulus was 



^ See http : / / www- xray . ast . cam . ac. uk/ papers/ dsdepro j / 

^ A fit was also done with the foreground column density al- 
lowed to vary, but the differences in the best-fitting parame- 
ters were negligible. We also note that the value of Ah de- 
termined by [Kalberla et al.| (2005) is ~ 10 per cent higher 
(Ah = 3.32 X 10^*^ cm~^). Fit results did not change significantly 
with Ah fixed to this value. 
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Figure 1. (a) Exposure-corrected, 0.5-3 keV image of HCG 62, smoothed with a gaussian of cr = (b) Residual image (in the sense 
of [model - image] / model) with contours, made using the exposure-corrected 0.5-3 keV image shown in panel a. The model was made 
using a Multi-Gaussian Expansion fCappe llari et a'L]|2006| >. The linear scale is identical to that of panel a. (c) Unsharpmask image of 
the core made by subtracting an image smoothed at 32^' from an image smoothed at (d) Unsharpmask image of the inner core made 
by subtracting an image smoothed at 5'' from an image smoothed at 



calculated as P = nkT, where we have assumed an ideal gas 
and n ^ 2ne. The entropy is defined as 5* = kTrie^^^. The 
cooling time was derived from the temperature, metallicity, 
and density using the cooling curves of Smith et al. ( ,2001 ), 
calculated using APEC. 

The resulting radial profiles of temperature, density, 
pressure, metallicity, and cooling time are shown in Figure |3] 
A number of features of these profiles are noteworthy. First, 
there appear to be two clear discontinuities in the temper- 
ature profile at r ^ 60'' and at r ^ 90^^ The relation of 
these discontinuities to the surface brightness edge seen at 
a similar radius in Figure [2] is discussed further in §3.1| 

The second notable feature in the radial plots is an off- 
centre peak (or, equivalent ly, a central dip) in the metallicity 
profile (this dip is also visible in the temperature map dis- 
cussed in {2.2.2). The metallicity has a value of ~ 0.3 solar 



in the innermost bin, rising to approximately solar metalli- 
cally at a radius of 20^^-30^^ then dropping to one- half solar 
at the group outskirts. Similar features have been noted in 
a number of groups and clusters (e.g., the Centaurus Clus- 
ter, ^^^^^^^^^ianjj2002). There has been some debate 
in the literature as to whether such features are real or not, 
and a number of possibilities have been proposed to explain 
them, such as resonance scattering (e.g., "Sanders & Fabian 
f2006) and modeling effects due to the presence of multi- 
phase gas (e.g.,|Buote Fabian|1998| |Buote|2000| [Molendi 



Gastaldeliol[2QQT[ 



Resonance scattering has been found to be generally 
insufficient to account for the observed metallicity deficits 
(Sanders k Fabian" '2006), and we do not investigate this 
possibility further. Instead, we examine the possibility that 
multi-phase gas in the core of HCG 62 has resulted in metal- 
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Figure 3. Radial profiles of temperature, density, pressure, metallicity, entropy, and cooling time (note that 10^^ equals 2.8 kpc). 



licit ies that are biased low. To this end, we followed the anal- 
ysis of [Kirkpatrick et al. (2009) and added in the inner re- 
gions a multi-phase cooling-flow component (MKCFLOW) 
with a low temperature fixed to 0.1 keV. The metallicity 
and high temperature parameters of the cooling-flow com- 
ponent were allowed to vary freely. Although an F-test does 
not indicate that the cooling-flow component improves the 
fit significantly, we do find that the metallically of the single- 
temperature component increases in the inner three regions 
with the addition of the cooling-flow component (see Ta- 
ble |2] and Figure [3| , largely erasing the central metallicity 
dip. This result is true for fits to both the projected and 
deprojected spectra, although the deprojected spectra show 
smaller metallicities in general (note, too, that the models 
generally fit much better to the deprojected spectra than to 



the projected ones). Therefore, it is possible that the central 
metallicity dip is due to our modeling of a multi-phase gas 
in the core with a simple single-temperature model. 



2.2.2 Temperature and Abundance Maps 

With the large number of counts available with the new 
data, it is useful to derive 2-dimensional maps of temper- 
ature and metallicity. To this end, bins were constructed 
for spectral extraction from the raw 0.5-7 keV counts im- 
age shown in Figure [l] using two binning algorithms: the 
Contbin algorithm of Sanders ( 2006 ) , which constructs bins 



that roughly follow contours in the counts image, and the 
WVT algorithm of |Diehl Statler] ( [2007t , which constructs 
bins using Voronoi tessellation. The bins were constructed to 
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Table 2. Parameters of the Spectral Fits in the Core. 



Annulus^ 


MetalUcity Xu 
IT IT+CF IT 


(dof) 

IT+CF 


Projected Spectra 


1 
2 

3 


0.49ig;5^ l-20ij-j^ 1.20 (672) 
0.60lg-Jg 1.63+°-^4 3.38 (670) 
0.95ig-;?o 2.42l2-;24 5 ^7 (g^g) 


1.09 (669) 
3.31 (667) 
5.04 (667) 


Deprojected Spectra 


1 
2 

3 


0.30+g-°J 0.93+°-^^ 1.45 (121) 
0.58ig-:°^ i.ii+o>i i_2g (121) 

^.79tlf2 l-68iJ:l^ 1-20 (121) 


1.29 (118) 
1.25 (118) 
1.13 (118) 



"See Figure |3] for inner and outer radii of the annuh. 



contain at least 2000 counts each (totaled over all three ob- 
servations) . Spectra and responses were extracted from each 
observation separately and were then fit jointly with an ab- 
sorbed single-temperature model in Sherpa. As before, the 
absorbing column density was fixed to the foreground hy- 
drogen column density of A^h = 3.03 x 10^° cm~^, and the 
redshift was fixed to z = 0.0137. 

The resulting maps of temperature and metallicity for 
each binning scheme are shown in Figure |4] Beyond the 
general trends seen in the radial temperature and metallicity 
profiles of higher temperatures and lower metallicities with 
increasing distance from the core, there is a clear asymmetry 
in both maps, with low-temperature, high-metallicity gas 
extending further to the south of the core than to the north. 
In the metallicity map, the asymmetry appears as a high- 
metallicity arc extending to the south and west of the core. 
This arc is discussed further in Section l4?2l 



2.3 Nuclear Point Source Emission 

NGC 4778, the BCG of HCG 62, is thought to host a low- 



luminosity AGN (L L AGN) at its core (Coziol et aL]|1998 
Shimada et al.|2000D . [Morita et aL] ( |2006 ) found no evidence 



for point-like X-ray emission at the location of the optical 
core of NGC 4778 (qj2ooo = 193.2738°, ^jsooo = -9.2040°) 



reported by [Zabludoff k Mulchaey (2000). The new data 



confirm this finding. However, Rafferty et al. (2008) report 
a somewhat different location for the optical core (q;j2ooo = 
193.2738°, (5j2ooo = -9.2040°), differin g hj ^ 2 arcsec from 
that reported by Morita et al. (2006). The reason for this 



difference is unclear. However, the Rafferty et al. (2008) lo- 
cation coincides precisely with possible point-like emission 
in the 2-7 keV band. To understand whether this emission 
could include a contribution from an AGN, we performed 
spectral fitting as follows. 

A spectrum was extracted from each observation in a 
circular region of radius 0.83^^ centred on the location of 
the optical core of NGC 4778 as reported by Rafferty et al.| 
(2008). The region was sized to enclose 90 per cent of the 



energy of an unresolved source at 1 keV at the position on 
the CCD of NGC 4778. As before, we extracted background 
spectra from the appropriate blank-sky background files in 
the same 0.83^^ region. 

The spectra from all three observations were fit simulta- 



neously between 0.5-7 keV in XSPEC with a model that in- 
cluded a thermal component (to account for IGM emission) 
and an absorbed power-law component (to account for AGN 
emission). As before, a WABS component was also used to 
account for Galactic absorption with the column density 
fixed to 3.03 x 10^° cm~^. A ZWABS component was used to 
account for absorption intrinsic to the AGN. The model was 
then defined as WABSx (MEKAL-hZWABSxPOW). Due to 
the low number of total counts (237, of which 4 are back- 
ground; in the 2-7 keV band, there are 22 source counts and 
1 background count), the Cash statistic was used on spectra 
binned to a minimum of 3 counts per bin (poor results were 
obtained when unbinned spectra were used). The power- law 
index was fixed to 1.8, typical of low- luminosity AGN (e.g., 
Ho et al.||2001||Terashima Wilson 2003). The metallicity 
of the MEKAL component was fixed to the value derived for 
the central bin from deprojection (0.3 Zq). The temperature 
and normalization of the MEKAL component were allowed 
to vary, as was the column density of the ZWABS compo- 
nent. The redshifts of the MEKAL and ZWABS components 
were fixed to ^ = 0.0137. 

Our fits confirm the presence of an X-ray AGN at the 
core of HCG 62. From the best-fitting normalization of the 
power-law component, we find an implied intrinsic X-ray lu- 
minosity for the AGN of Lo.5-7 keV = (1.5ti:o)xl0^^ ergs"^ 
(3-cr errors are quoted). The intrinsic column density is only 
poorly constrained, with the best-fitting value being A^h = 
0.6to 5 X 10^^ cm"^ (1-cr errors) and A^h < 6.0 x 10^^ cm"^ 
(3-cr upper limit). Therefore, the spectrum shows evidence 
for moderate absorption, as seen in a number of LLAGN 
(e.g. 



Terashima & Wilson||2003 ). 



Using the Ha flux implied by the measurements of |Co-| 
ziol et al. ( |1998| and Shimada et al. ](|2000|) for the nucleus of 
NGC 4778 (Fhc = 7.3 x 10~'^ erg cm"^ s"^), we can com- 
pare the ratio of the X-ray luminosity to the Ha luminosity 
to that of typical LLAGN. Using the 2-10 keV luminosity 
implied by our fit of L2-10 keV 
find a ratio of L2-10 kev/-t/Ha 
typical ratios in LLAGN of ^ 
2000||Ho et aL||200T 1). 



(I.ld=0.4)xl0^^ ergs-\ we 
' 2.7 i: 1.0, consistent with 



1-2 (e.g., j Terashima et al. 



Lastly, we find evidence for variability of the 2-7 keV 
nuclear X-ray emission over the period spanned by our ob- 
servations 9 yr). We find background-subtracted count 
rates in the 2-7 keV band of (4.4 ± 3.1) x 10~^ cuts s"^ 
for ObsI D 921 (taken in early 2000; errors are calculated 



following Gehrels 1986) and (15 ± 5) x 10"^ cuts s"^ and 



(21±7)xl0"^ cuts s"' for ObsIDs 10462 and 10874, respec- 
tively (both taken in early 2009). Therefore, it appears the 
AGN has brightened by a factor of ~ 3 during this 9-year 
period. 



3 THE AGN OUTBURST 

As previous studies of HCG 62 have reported, the AGN has 
interacted strongly with the IGM in the core of the group. 
The large outer cavities visible in Figure ^ are the most 
prominent feature of this interaction, but there is also the 
inner surface brightness edge noted in §2.1| that, given its 
proximity to the cavities, is likely associated with them. 

In this section, we examine the structure, energetics, 
and content of the cavity system. We first examine the sur- 
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Figure 4. Left: Temperature maps made using contour binning {top) and using WVT binning [bottom) with contours from the residual 
image (Figure [l]) overlaid. Errors on the values in each bin are typically 2%-4%. Note the presence of cooler gas to the inside of the 
surface-brightness edge (indicated by the southern-most contour) and warmer gas to the outside, a typical signature of a cold front. 
Right: Metallicity maps (relative to solar) with contours from the residual image overlaid. Errors are typically 20%-30%. 



face brightness edge that encompasses the outer cavities to 
determine whether this feature could be due to a weak shock 
(which would contribute to the total energetics of the out- 
burst) or might instead be due to rims of cool material sur- 



rounding the cavities such as those as seen in A2052 (Blan- 



ton et al.||2QQ3| ). We then estimate the total outburst ener- 
getics and place limits on the contents of the cavities. 



3.1 The Surface Brightness Edge 

To determine whether the surface brightness edge surround- 
ing the cavities (see Figure |5| is consistent with a shock 
discontinuity, we fit a broken power-law model to the az- 
imuthally averaged surface brightness profile shown in Fig- 
ure [6] In the fit, we chose the power-law indices inside and 
outside of the edge to match our observed surface brightness 
profile in those regions, and solved for the best-fitting radius 
and density jumpy for the discontinuity. As the power-law 
density profile is only an approximation to the true one, we 



fit over a limited range of radii near the location of the edge. 
Lastly, a background of 7.9 x 10~^° counts s~^ cm~^ pixeP^, 
determined by fitting a double-/^ model plus a constant back- 
ground to the surface brightness profile (see Table [l]), was 
subtracted before fitting was performed. 

The resulting best-fitting discontinuity parameters are 
given in Table [S] and the fit is shown in Figure |6] This dis- 
continuity model fits acceptably (given the presence of sig- 
nificant non-azimuthally symmetric structure, such as the 
cavities, that distorts the surface brightness profile some- 
what) and implies a density jump of 1.45 across the edge 
and a Mach number of 1.33. 

However, if we fit to surface brightness profiles derived 
in subsectors (shown in Figure [5|, we do not find consistent 
fit parameters (see Table p| Fl Although sectors 3 and 5 are 



^ We also find this behavior if elliptical regions, with ellipticity 
and position angle matched closely to those of the edge (e = 0.15 
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Table 3. Parameters of the Broken Power-law Models. 



Region" 



Break Radius 
(pixels) (kpc) 



Density Jump 



Inner ec 


ge, all sectors 


127.5 


17.5 


1.45 ± 0.04 


2.39 


Inner ec 


ge, sector 1 


128.4 


17.7 


1.82 ± 0.17 


1.32 


Inner ec 


ge, sector 2 


133.2 


18.3 


1.93 ± 0.18 


2.12 


Inner ec 


ge, sector 3 


139.3 


19.2 


1.98 ± 0.18 


9.28 


Inner ec 


ge, sector 4 


103.9 


14.3 


1 9-1+0.15 


2.39 


Inner ec 


ge, sector 5 


103.8 


14.3 


1.50 ± 0.10 


3.89 


Inner ec 


ge, sector 6 


210.4 


29.0 


1.07 ± 0.4 


2.01 


Outer ec 


ige, sector 5 


256.6 


35.4 


1.48 ± 0.17 


1.11 



"See Figure [S] for region definitions. 

clearly affected by the cavities and fits to the profiles in 
these sectors may not be meaningful, fits to the other sectors 
should be consistent with one another. 

Instead, there is a large range in the best-fitting radii 
and density jumps among these sectors, with those in sec- 
tors 4 and 6 differing the most. Due to symmetry, one might 
expect that the fits in sectors 1 and 4, located on opposite 
sides of the core, would be similar; however, the best-fitting 
density jumps for these sectors differ at the 3-cr level. Fur- 
thermore, in sector 6 there does not appear to be any signif- 
icant edge in the surface brightness profile at the expected 
radius. 

Although substructure in the core makes it difficult to 
obtain reliable fits, our analysis implies that a single discon- 
tinuity model cannot readily account for the surface bright- 
ness edge around the entire circumference. Furthermore, the 
inner temperature jump seen in Figure |3] corresponds closely 
to this edge and indicates that the temperature behind the 
edge is significantly lower than that beyond the edge. These 
results disfavor the shock interpretation. An alternative ex- 
planation for this edge is that it is due to the presence of 
rims associated with the cavities. We examine this possibil- 
ity further in §3.2[ 



3.2 Energetics 

Measurements of X-ray cavities have for the first time al- 
lowed a strong lower limit to be placed on the total jet power 
of AGNs. Figure [l]: shows an unsharp mask of the central 
region of HCG 62, showing the large cavities to the north 
and south of the core. These cavities are filled with radio 
emission ( Gitti et al.||20ld ). There is also evidence of small 
inner cavities, oriented roughly north-south. The jet power 
of the outburst may be estimated from the enthalpies and 
buoyant rise times of the cavities. 

To quantify the significance of the cavities, we can com- 
pare counts in the cavity regions shown in Figure [l] relative 
to those in undisturbed regions of the same size and at the 
same radius from the core. We find that the outer cavities 
represent deficits of ~ 28 per cent (2062 it 117 counts) for the 
southern cavity and ~ 14 per cent (1513 ± 142 counts) for 
the northern cavity (this value includes the spur of bright 
emission in the northern cavity visible in Figure [l]:). For the 



and 6 = 28° east from north), are used in the derivation of the 
surface brightness profile. 




Figure 5. The exposure-corrected 0.5-3 keV counts image with 
the six sectors used for fitting of the shock models to the the 
inner surface brightness edge. The inner and outer boundaries 
of the sectors correspond to 50 and 200 pixels (25^^ and 98^^), 
respectively, the approximate range over which fitting of the inner 
edge was done (fitting of the outer edge was done in sector 5 
between 120 and 500 pixels). 



inner cavities, visible in Figure [T]l, the deficits are 16 per 
cent (121 ± 37 counts) for the southern cavity and ^17 per 
cent (198 lb 39 counts) for the northern cavity. Therefore, the 
significance of the inner cavities is much lower than that of 
the outer cavities, and, consequently, we examine the large 
outer cavities in much greater detail. 

To make the best possible estimate of the energetics 
of the outer cavities, we used simulated images to recover 
the cavity sizes and locations. The large-scale cluster emis- 
sion was modelled with the 2-D double /3-model described 
in §2.1| the cavities as empty spheres, and the rims as uni- 
form shells of emission that fully enclose the cavities. The 
0.5-3 keV exposure-corrected image was binned using the 
WVT binning method of |Diehl St atler f2006^ to a S/N 
of 25 per bin. The same binning scheme was then applied 
to the model image at each iteration of the fitting process, 
which consisted of a least-squares fit performed with the 
TNMIN function of MPFIT (Mar kwardt|2009t . The spur of 
bright emission in the northern cavity was masked, as was 
the bright emission that runs between the cavities, perpen- 
dicular to the cavity axis (see Figure [l]:). A summary of the 
best-fitting cavity parameters is given in Table ^ and the 
residual image is shown in Figure [7| For this fit, we obtained 
= 3.2 (this includes contributions from non-azimuthally 
symmetric regions of the image, such as the arc of excess 
emission, that were not well fit by the 2-D double /3-model). 

It is clear from Figure |7| that the cavity regions and in- 
ner surface brightness edges are well fit by the model. How- 
ever, there remain residuals in the cavity regions that may 
indicate that they are distorted somewhat from spheres. In 
the northern cavity, some of the residuals may be due to 
the spur of bright emission noted above. The nature of this 
spur is unclear. One possibility is that it is cool material up- 
lifted by the cavity in its wake, and, in this case, the bubble 
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Figure 6. Top: The azimuthally averaged, 0.5-3 keV surface 
brightness profile across the inner surface brightness edge (be- 
tween radii of 50 and 200 pixels or 6.9 and 27.5 kpc). The best- 
fitting broken power-law model is shown as the solid line. Bottom: 
The 0.5-3 keV profiles and best-fitting broken power-law models 
across the inner surface brightness edge for the six sectors defined 
in Figurejs] A background of 7.9 x lO"-*^*-* counts s"-*^ cm~^ pixel"^ 
(determined from the double-/? model fit described in ^2.1| was 
subtracted from both plots. 



is probably nonspherical. However, deriving the true cavity 
geometry would be difficult due to projection effects and is 
beyond the scope of this work. 

Interestingly, the total emissivity in the rims of the best- 
fitting model is very close to that displaced by the cavities 
(within 10 per cent). This similarity could imply that the 
rims are composed of this displaced material. However, it 
could also be simple coincidence. A similar analysis for a 
sample of systems with cavities would be useful in under- 
standing whether such a situation is common or not. We 
note that the inclusion of rims does not change the best- 
fitting cavity properties appreciably, but it does largely ac- 
count for the inner surface brightness edges discussed in §3.1| 
However, given the complex structure in the core, it is dif- 
ficult to say conclusively whether the rim interpretation or 
the weak shock interpretation is the correct one (some com- 
bination of structure from both a rim and a shock is also 
possible). 

The best-fitting values of the cavity sizes and locations 
give a total cavity enthalpy {H = ApV^ where p is the average 



Table 4. Parameters of the Best-fitting Cavity Model. 

r R (f)'' 0^ Rim width 

Cavity (kpc) (kpc) (deg) (deg) (kpc) 



North 
South 



5.5 
4.3 



9.3 
6.9 



90 
90 



135 
288 



5.5 
4.4 



0.89 
0.95 



"Angle to the line of sight of a line connecting the cavity centre 
to the group core. 

^Angle in the plane of the sky, increasing counter clockwise north 
from west, for a line connecting the cavity centre to the group 
core. 

^Ratio of the total emissivity in the rims relative to the emissivity 
displaced by the cavities (i.e., a ratio of unity would imply equal 
emissivity in the rims and in the region displaced by the cavity). 

pressure at the radius of the center of the cavity, calculated 
from the pressure profile shown in Figure |3] V is the vol- 
ume, and assuming a relativistic plasma fills the cavities) 
for the outer cavities of iftot ~ 3.6 x 10^^ erg (see Table pi 
for the derived cavity properties). This value is somewhat 
larger (by a factor of two) than that found by other stud- 
ies (e .g., |B!rzan et aL]|2004| [Morita et aT1|2006| |Gitti et al. 



2010| ), although it is a factor of two smaller than the value 



derived by Gitti et al. (2010) using cavity sizes estimated 
from the radio emission. For the small inner cavities, we es- 
timate (by fitting ellipses to them by eye) a total enthalpy 
of ~ 1.1 X 10^^ erg. Using the buoyant rise times as the ages 



2004), we find 



(tbuoy) of the cavities (e.g., iBirzan et al 
total mechanical power (Hi 

Lmech 8.6 X 10^^ erg s , ^ 

luminosity of 1.8 x 10^^ erg s"^ ( |Rafferty et al.||2006 ) 



tot/tbuoy) for all four cavities of 
a factor of ^ 5 above the c ooling 



3.3 Cavity Contents 



erg s 



(Gitti 



Given the total radio luminosity of 3 x 10^^ 
et al. |2010 ), the implied radiative efficiency (defined as 
^rad/[^rad + ^mech]) Is ^ 4 X 10~^. Under the assumptiou 
of equipartition, the requirement of pressure support con- 
strains the content of the lobes. For the observed radio lu- 
minosity, the relativistic electrons in the lobes can provide 
only a tiny fraction (< 0.2 per cent) of the particle pressure 
requi red ( [Morita et al]|2006| [Birzan et aLl|2008| [Gitti et al 
20101. 



This result implies a large population of heavy particles 
(such as cosmic rays) or hot thermal gas in the lobes. For 
a given cavity geometry and background IGM distribution, 
one can use the surface-brightness residuals in the cavity re- 
gions to place limits on any thermal gas in the cavities (e.g,. 
Schmidt et al.[[2002| [Morita et al.[[2006[ ). However, the true 
cavity geometry is difficult to constrain (as are deviations 
in the gas distribution from the beta model atmospheres 
assumed in most analyses), and arguments which rely criti- 
cally on knowledge of the true cavity geometry and gas dis- 
tribution may not provide useful constraints. For example, 
if the cavities are elongated along the line of sight (relative 
to spherical cavities) , a thermal gas inside the cavities could 
be present but result in little to no residual emission when a 
spherical cavity model is used. To minimize these issues, we 
adopt a method similar to that used by [Sanders &: Fabian[ 
([2007,) that uses spectral constraints and depends much less 
on knowledge of the true cavity geometry. 

To this end, we extracted spectra from circular regions 
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Table 5. Cavity Energetics. 



Cavity 



(10 erg cm ^) 



V 

(10^^ cm3) 



(1056 grg) 



^buoy 

(10^ yr) 



(10^2 erg s-i) 



Inner N 
Inner S 
Outer N 
Outer S 



9.7± 1.1 
9.7± 1.1 
2.9 ±0.2 
3.3 ±0.2 



0.12 
0.16 
20.8 
9.5 



0.48 
0.6 
24.0 
12.4 



0.31 
0.30 
1.7 
1.3 



0.49 
0.63 
4.5 
3.0 




Figure 7. Residual images (in the sense of [image - model] /image) of the core, made using the 0.5-3 keV image binned to a S/N of 25 
per bin, for a model that includes only a 2-dimensional double-Beta component {left) and for the best-fitting model that also includes 
empty spherical cavities surrounded by uniform spherical shells {right). Masked regions are shown in black. The outer cavities and inner 
surface brightness edges are marked. 



of radius in the deepest regions of the outer cavities, 
where any contribution from thermal gas inside the cavities 
relative to gas outside the cavities should be greatest. The 
spectra, responses, and backgrounds were made as described 
in §2.2| We fit the spectra with a model that includes two 
single-temperature components. The first component, with 
the temperature fixed to the average of that at the radius 
of the cavity, accounts for the projected emission. The sec- 
ond component accounts for any thermal emission inside the 
cavity, and its temperature is fixed to a value that was var- 
ied from 1.5 keV (the maximum temperature seen in the 
temperature profile of Figure [s]) to 5 keV. The metallicities 
were fixed to the the value obtained from deprojection at 
the cavity radius. Lastly, a WABS component was added to 
account for Galactic absorption, with the A^h column den- 
sity fixed to 3.03 x 10^*^ cm~^, the value used throughout 
this paper. The full model is then WABS x (MEKALi + 
MEKAL2). 

Assuming that any emission from gas at temperatures 
above 1.5 keV comes from gas inside the cavity that is in 
pressure balance with the gas surrounding the cavity, we 
can place constraints on the volume-filling factor of this 
gas using upper limits on the normalization of the second 
MEKAL component. For the temperature and density of 



the surrounding gas, we use the deprojected values (from 
the deprojection analysis discussed in ^2.2.1) at the radius 



of the cavities. We assume that the cavities in the regions 
from which spectra were extracted are cylinders of length 
11 kpc and 8.5 kpc for the northern and southern cavities, 
respectively (corresponding to the diameters recovered from 
our fits described in ^3.2). We find that a volume- filling gas 



with a temperature less than 4.3 keV is excluded by our 
spectral fits in the southern cavity (see Figure |8| . Fits to 
the northern cavity did not provide useful constraints (i.e., 
we could not exclude the presence of volume-filling gas of 
any temperature above 1.5 keV), probably due to the lower 
contrast any such gas would have at the larger projected 
radius of this cavity. We note that differences of 20 per cent 
in our adopted cavity diameter result in changes of 10 per 
cent in the temperature limit (e.g., a diameter 20 per cent 
smaller results in a limiting temperature of ~ 3.9 keV). 

Our limit of 4.3 keV on the temperature of volume- 
filling gas in the southern cavity is somewhat higher than 



the value of 3 keV determined by Morita et al. ( 2006 ) under 



the assumption of spherical cavities. However, it does not 
of course rule out the possibility that significant pressure in 
the cavities could be supplied by a gas with a temperature 
above this limit. Although deeper observations would allow 
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Figure 8. The limit on the volume-fining fraction versus the 
temperature of any thermal gas in the southern cavity. 



this temperature limit to be increased, it is impossible to 
rule out gas at arbitrarily high temperatures (due to the 
scaling of emissivity with temperature for a gas in pressure 
balance with its surroundings). 



4 A RECENT MERGER 

In this section we examine evidence for a recent merger 
in the surface brightness distribution and temperature and 
metallicity structure of the IGM. In particular, we investi- 
gate whether the outer edge noted in §2.1| is consistent with 
a weak shock (possibly generated by an AGN outburst) or 
rather with a cold front (generated by a merger), and we 
examine how the spiral-like excess of emission visible in Fig- 
ure [TJ) relates to this edge and to the merger scenario. 



4.1 The Outer Edge: Shock or Cold Front? 



As noted by |Gitti et al. ( 2010 ) , there is a clear surface bright- 
ness edge ~ 2 arcmin (~ 30 kpc) to the south of the core. It 
was posited by Gitti et ah] ( |2010| ) that this edge is associated 
with a weak shock, such as the one seen in MS0735+7621 



( McNamara et al.^20 05), that was possibly generated by the 
AGN outburst which created the cavities and radio lobes. 
In this section, we investigate whether this edge is consis- 
tent with a discontinuity formed by a weak shock and we 
examine the temperature structure of the gas near the edge. 

We first fit the surface brightness profile across the edge 
with the broken power-law model discussed in §3.1| We ex- 
tracted the profile from the region of the outer edge only 
(sector 5, shown in Figure|5]), and, as before, we subtracted a 
constant background of 7.9 x 10~^° counts s~^ cm~^ pixel"^. 

We find a good fit to the observed profile (see Figure|9|, 
with discontinuity parameters (see Table [sj that match well 
with those found by Gitti et ah] ( |2010| (although our best- 
fitting model implies a density jump 10 per cent lower). 
Therefore, the edge in this region is consistent with the dis- 
continuity expected to accompany a weak shock. 

However, the temperature structure of the IGM in this 
region does not show the expected signature of a shock: 
namely, the temperature immediately behind the shock 
should be higher than that immediately in front. Instead, 
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Figure 9. Surface brightness profile (0.5-3 keV) and best-fitting 
shock model for the outer surface brightness edge (note that 
100 pixels equals 13.8 kpc). A background of 7.9 x 10"-*^° counts 
cm~^ pixel"^ (determined from the double-/? model fit de- 
scribed in §2.1|) was subtracted from the plot. 



the temperature maps (see Figure |4| show that the temper- 
ature outside of the edge appears to be higher than that 
inside. Therefore, the temperature structure suggests that 
the discontinuity may be a cold front, and we examine this 
possibility further in the following section. 



4.2 



The Spiral-like Excess and High Abundance 
Arc 



The arc of high metallicity gas not ed in §2. 2. 2 1 was d iscovered 
in the original Chandra data by Gu et al.| (|2007 ) and was 
confirmed by a reanalysis of those data by Gitti et al. ( 2010 ). 



In the new data, this arc is seen to correspond closely to a 
region of excess surface brightness (relative to the model of 
the large-scale emission described in ^2.1 see Figure [I). The 
southern edge of this arc defines the outer surface brightness 
edge discussed in §4.1| and hence the arc and the edge ap- 
pear to be related. Although a weak-shock model fits the 
edge well (Gitti et al. 2010| ), it is difficult to understand 



the temperature and metallicity structure of the arc in this 
context. 

Alternatively, the arc and its edge could be indications 
of a recent merger, as suggested by Gu et al. (2007|. In- 
deed, these structures appear to match well with those seen 
in simulations of mergers in cooling flows (e.g., Ascasibar 
Markevitch 2006). In this scenario, the merger induces 



sloshing in the cool core, uplifting material from the cen- 
tre, and creating a low-temperature (and presumably high- 
metallicity) spiral-like feature and associated cold front (the 
outer edge). Additional support for the merger scenario in 
HCG 62 is provided by the kinematic study of Spavone e t al.| 
(|2006|), who found evidence that NGC 4778 has undergone 



a recent interaction, possibly with NGC 4761. 

Therefore, although it is not possible to distinguish ab- 
solutely between the weak shock scenario and the merger 
scenario as the cause of the arc and surface brightness edge, 
the evidence obtained from the new Chandra data appears 
to favor the merger scenario. 



© 0000 RAS, MNRAS 000, 000-000 



12 Rafferty et al. 



5 DISCUSSION 

5.1 The Recent Merger 

We have presented evidence for a recent merger in the 
structure of the IGM. Our results suggest that the high- 
metalhcity arc to the south of the core could be gas uplifted 
by sloshing induced by the merger. Alternatively, the high- 
metallicity arc material could have been stripped from the 



merging galaxy, as postulated by Gu et al. (2007). Lastly, 



some of the metals in the arc could have been uplifted not 
by sloshing induced by the merger, but instead by the AGN, 



as seen in many clusters (Kirkpatrick et al. 112011) 



However, the values of metallicity in the arc appear to 
exceed the central metallicity. It would therefore be diffi- 
cult to understand how this gas could originate in the core. 
However, as we showed in §2.2| the true central metallicity 
could be higher than a single-temperature model indicates 
due to the presence of multi-temperature gas in the core. 
Indeed, the inclusion of a multi-temperature component (an 
MKCFLOW component) to fits in the core does result in 
central metallicities that are consistent with those in the 
arc. Therefore, in this scenario, the material in the arc could 
have been uplifted from the core, either by sloshing or by 
the AGN. 



The simulations of Ascasibar & Markevitch ( 2006 ) show 



that sloshing is capable of uplifting cool material from the 
core and can result in structure similar to that seen in the 
surface brightness and temperature maps of HCG 62 (i.e., 
cold fronts and spiral-like arcs of lower-temperature mate- 
rial). However, it is not clear that AGN activity could easily 
replicate this structure, as such activity generally appears 
to distribute the metals along the cavity or radio jet axis 



in a fairly symmetrical way (e.g., Kirkpatrick et al.||201l| . 
Furthermore, there is no evidence of enhanced metallicity in 
the direction of the northern lobe and cavity. Therefore, it 
appears unlikely that AGN activity alone could account for 
the observed metallicity distribution. 

Whether the metals originate in the group core or were 
instead stripped from the merging galaxy is more difficult 
to constrain. However, we note that the majority of ma- 
terial in the group with the highest metallicities would be 
expected to occur in the core where it would undergo the 
most enrichment. Additionally, there is no evidence of high- 
metallicity gas at the location of the probable merger can- 
didate, NGC 4761, suggesting that the stripping must have 
been very efficient. It therefore seems probable that sloshing 
is the primary mechanism by which the high-metallicity arc 
was formed. 

Such sloshing could have a significant heating effect on 
the cool core, as cool gas is mixed with hotter gas and turbu- 
lence increases (e.g., |Sharma et al.|2009||ZuHone et al.|2010| ). 
However, in contrast to the cavities, any heating from such 
mechanisms is difficult to quantify. Furthermore, increased 
turbulence could conceivably affect entrainment of gas by 
the radio lobes, which might explain in part the high in- 
ferred values of the ratio {k) of the total particle energy to 
the energy in electrons for the outer cavities of HCG 62 (on 
the order of hundreds; [Birzan et al.|2008| [Gitti et al.|2010t . 
In this scenario, one might expect that the inner cavities 
would have much lower values of /c, as they would presum- 
ably have had less time to suffer the effects of entrainment. 



High resolution radio images of the inner cavities are needed 
to test this possibility. 



5.2 The AGN Outburst 

The cavity system in HCG 62 has been well studied, but 
there remain some poorly understood aspects. For instance, 
the deep radio image at 235 MHz presented in [Gitti et al.| 
(2010) clearly shows extensions of emission beyond the 



boundaries of the X-ray cavities. These extensions are spec- 
trally steep, as they are not detected in the deep 610 MHz 



image of Gitti et al. (2010). We have examined the regions 



where these extensions lie and find no evidence for signif- 
icant X-ray structure that might be associated with them 
(such as older "ghost" cavities) . It is possible that such emis- 
sion leaked from the cavities, as seen for example in M84 



(Finoguenov et al. 2008). However, given the small size of 



the extensions and their large distance from the core, the 
deficits produced by cavities would likely be undetectable in 
the current image. 

If leakage of radio plasma out of the cavities did occur, 
leakage of thermal gas into the cavities might then be ex- 
pected. Again, such a scenario is consistent with the high 
values of k inferred for the cavities. If such gas has leaked 
into the cavities, it presumably has either been heated to 
temperatures above the limit of ~ 4.3 keV (for the south- 
ern cavity) or fills only a fraction of the cavity volume (e.g., 
^ 20 per cent for a gas at 2 keV). 

As for the energetics of the AGN outburst in HCG 62, 
we have shown that only ^ 20 per cent of the total cav- 
ity enthalpy is sufficient to balance cooling if it is injected 
on time scales similar to the buoyancy time scales inferred 
for the outer cavities. The remaining 80 per cent presumably 
goes (eventually) into net heating of the IGM. The net heat- 



ing from the cavities is then at least 0.8i/t( 



3 X 10^' erg. 



Given the total gas mass of 1.6 x 10 Mq within a radius of 
300 kpc ( [Morita et al.|2006| ) , the current outburst is supply- 
ing O.SHtot/J^rrip/Mgas ~ 6 x 10~^ keV per particle of heat 
(where nip is the proton mass and we have adopted a mean 
molecular weight, /x, of 0.59). Assuming a typical inter-burst 
timescale of ~ 1.5 x 10^ yr and summing this heat over the 
life of the group (^ 7.7 x 10^ yr, the lookback time to z = 1), 
we find a total heating by the AGN of ~ 0.3 keV per particle. 
This value is somewhat higher than the average value of 



0.17 keV per particle inferred by Ma et al. (2011) from the 
radio luminosities of a sample of 169 clusters. It is signifi- 
cantly lower than the heat required to break self-similarity 
1-3 keV per particle; Wu et al.|2000 ). However, the heat- 
ing rate could well have been higher in the past, when the 
supermassive black hole at the core of the group was likely 
experiencing the bulk of its growth. Indeed, studies that es- 
timate the average heating indirectly, such as from estimates 
of the black hole mass (e.g., Wang et al.||2010 ) or from the 
gas profiles (e.g., Mathews & Guo^^n|, find much larger 
values of heating, implying higher heating rates at earlier 
epochs. 



6 CONCLUSIONS 

We present results of an analysis of new Chandra data of 
HCG 62. Our main findings are as follows: 
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(i) We refine previous measurements of the AGN cavities 
and confirm the extremely low radiative efficiency of the 
outburst. This low efficiency constrains the pressure support 
from electrons in the cavities to be much less than needed 
to balance the external pressure of the IGM, implying that 
the cavities contain either thermal gas (possibly entrained 
by the lobes) or a signiffcant population of cosmic rays. 

(ii) Using spectral fits, we constrain the temperature of 
volume-filling thermal gas in the southern cavity to be kT > 
4.3 keV. However, gas at lower temperatures with a lower 
volume filling fraction is allowed by our fits. 

(iii) For the first time, we detect the presence of an X-ray 
AGN at the core of HCG 62. We constrain the X-ray lumi- 
nosity of the central AGN to be Lo.5-7 kev = (^-^-I'.o) x 
10^^ erg s~^ (3-cr errors). The X-ray spectrum shows evi- 
dence of moderate absorption (A^h ~ 0.6 x 10^^ cm~^). 

(iv) We find evidence of a recent merger in the surface 
brightness, temperature, and metallicity structure of the hot 
gas. Such a merger may have contributed to the heating of 
the cool gas in the core, although the AGN cavities alone 
represent more than enough power to offset cooling. 

In the future, high resolution radio observations would be 
useful in verifying the inner cavities detected in the new data 
at low significance and in constraining their properties. 
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